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Summary 


This  report  summarizes  pi  ogress  on  three  currently  active  efforts 
in  the  investigation  of  laser  beam  scintillations  due  to  atmospheric  turbulence. 
These  topics  are  (1)  multiwavelength  scintillation  over  a  long  horizontal  path 
with  a  very-high  integ~ated-path  turbulence  level;  (2)  finite-beam  or  trans¬ 
mitter-aperture  effects  including  beam  wander,  spread,  and  scintillation;  and 
(3)  turbulence  intermittency  effects. 

The  long  path  results  are  presented  in  detail.  Saturation  of  scintil¬ 
lations  was  observed  at  10.  6  pm,  and  the  behavior  far  into  saturation  was 
explored  at  4880  &.  The  evolution  of  very  small  and  very  large  scintillation 
scale  sizes  was  observed,  as  predicted  in  a  recent  theoretical  treatment.  The 
effects  of  atmospherically-induced  beam  wander  and  spread  on  mean  target- 
illumination  are  explored  theoretically,  as  are  the  fading  effects  owing  to 
combined  wander  and  scintillations.  Certain  inconsistencies  in  the  literature 
pertaining  to  these  effects  are  resolved,  and  the  advantages  to  be  expected 
from  an  ongoing  beam-wander -cancellation  experiment  are  described.  Turbu¬ 
lence  intermittence  is  modelled  in  terms  of  discrete  slabs  of  turbulence,  and 
inner  scale  effects  are  analyzed.  The  intermittency  problem  is  redefined  and 
experimental  plans  are  described. 
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I.  Introduction 


During  the  reporting  period,  work  was  eonducted  on  long-path 
scintillations,  finite-beam  effects,  and  turbulence  intermittent  effects. 

and  th  ThC  l0ng'Path  eXPeHmen‘S  have  completed,  the  data  analysed, 
-d  the  results  interpreted.  This  effort  is  described  in  detail  in  Section  II 

ery  interesting  new  results  were  obtained  regarding  saturation  at  lOum  ' 
wavelengths,  the  evolution  of  new  scintillation  scale  lengths  at  large  inte 
grated -path  turbulence  levels,  and  the  behavior  far  into  saturation 

to  section  III,  the  theoretical  predictions  of  finite-beam  effects 
are  reviewed  in  detail,  including  beam  wander,  instantaneous  spread,  and 
scintillation.  Apparent  inconsistencies  and  contradictions  in  various  treat¬ 
ments  in  the  literature  are  resolved,  and  die  results  are  applied  to  the 
specification  of  important  new  experiments  in  this  area.  These  include  die 

cancellation  of  beam  wander  through  the  use  of  a  reciprocal  property  of  tur- 
bulence  propagation. 

with  r  dI”tSeCli°n  IV’  redef“e  tUrbUle"Ce  totermittency  problem 

6  Potion  effects,  and  present  some  pertinent  analytical 

results.  y 
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II.  Long -Path  Scintillations 

The  purpose  of  these  experiments  was  to  measure  the  properties  o£ 
scintillations  over  the  longest,  lowest  path  which  can  be  reasonably  achieved, 
in  order  to  encounter  the  highest  integrated-path  turbulence  level  possible. 
The  experiments  were  conducted  with  simultaneous,  coincident  virtual-point- 
sources  at  visible  and  middle -infrared  wavelengths.  The  measurements  are 
itemized  in  Table  I,  and  include  turbulence  strength,  meteorological  param¬ 
eters,  and  scintillation  statistics.  The  turbulence  and  scintillation  quantities 
are  defined  in  Ref.  1,  and  the  general  instrumentation  and  field  facilitiei  are 
described  in  Refs.  2  and  3,  respectively. 

The  experimental  parameters  are  summarized  in  Table  II.  Due  to 
the  fact  that  the  experiments  took  place  over  such  a  long,  uniform  path  near 
the  ground,  long-term  vertical  beam -bending  due  to  thermal  gradients  was 
significant.  Thie  was  manifested  as  a  mirage  effect  during  early  morning 
hours,  and  conversely  as  a  looming  horizon  during  afternoon  or  high-turbu¬ 
lence  hours.  As  a  result  of  the  latter,  it  was  necessary  to  utilize  elevated 
transmittere  and  receivers  to  maintain  an  unobstructed  path  during  peak 
turbulence  periods.  Significant  difficulties  were  experienced  with  laser  relia¬ 
bility,  such  that  reliable  data  were  not  taken  until  the  final  few  high-turbulence 

days  of  summer.  However,  the  quantity  of  data  is  considered  sufficient  for 
the  conclusions  drawn  below. 

Major  objectives  of  these  experiments  were  the  measurement  of 
saturation  of  log  amplitude  variance  and  related  effects  at  the  10.  6  urn  wave¬ 
length.  the  investigation  of  scintillation  statistics  at  a  visible  wavelength  far 
into  the  saturation  region,  and  the  examination  of  large-integrated-path  tur¬ 
bulence  effects  on  scintillation  correlation  scales  and  hence  receiver  aperture¬ 
smoothing.  In  the  following  sections,  we  discuss  log  amplitude  variance, 
covariance,  scintillation  spectra,  probability  distributions,  and  receiver 
aperture  smoothing  respectively. 


TABLE  I.  Experimental  Measurements. 

Strength  of  turbulence  (C^  )  from  microthermal  probes 

Vertical  temperature  gradient 

Wind  velocity 

Log  amplitude  variance 

Log  amplitude  covariance 

Log  amplitude  probability  distribution 

Scintillation  spectrum 

Receiver  aperture  smoothing 

TABLE  II.  Experimental  Parameters. 

Path  length:  6.  0  km 

Path  description:  Farmland,  flat  to  within  +  0.  5  m 
Transmitter  height:  6.  1  m 
Receiver  height:  3.  5  m 

Wavelengths  (simultaneous,  coincident):  4880  X,  10.  6  um 

Transmitter  beam  configuration:  Virtual  point  sources 
(Fresnel  number  <10“3) 

Receiver  aperture:  6  mm 

Receiver  dynamic  range:  >80  dB 

Receiver  averaging  time:  60  sec 

Aperture -smoothing  receiver:  0.6-32  cm 

Resolution  of  spectral  measurements :  1  Hz 

Receiver  band>yidth:  1  kHz 

Microthermal  probe  separation  for  C  2;  10  cm 

n 

Microthermal  probe  height:  1.8  m 
Microthermal  averaging  time:  300  sec 


A.  Log  Amplitude  Variance 

2 

Experimental  log  amplitude  variances  ( cr  )  were  obtained  from 

E 

probability  distributions.  A  typical  diurnal  behavior  for  a  cloudless  day  is 

shown  in  Figure  1.  The  unusually  high  values  of  c r^,  ^  for  4880  X  in  the  early 

2 

morning  were  due  to  poorly  developed  turbulence  structure,  and  there  is 

evidence  that  the  presence  of  light  ground  fog  during  that  period  further 

increased  these  anomalous  fluctuations.  The  data  show  that  for  this  very 

long,  low  path  configuration  the  variance  at  4880  X  was  nearly  always  saturated. 

As  discussed  below,  the  lower  extreme  of  these  variances  is  not  indicative  of 

an  asymptote  for  high  turbulence.  Outside  of  the  anomalous -turbulence 

period,  the  10.  6  qm  variances  were  seen  to  saturate  in  mid-morning. 

The  variances  measured  at  10.  6  (a,m  are  shown  in  Figure  2  vs. 

the  strength  of  turbulence  taken  at  a  height  of  1.  8  m.  For  the  first  time, 

saturation  is  clearly  indicated  at  this  wavelength.  As  discussed  below,  the 

abscissa  is  not  highly  meaningful  without  a  correction  due  to  variable  beam- 

refraction,  The  saturation  level  appears  to  be  somewhat  below  the  0.  6  value 

1  2 

typical  of  shorter  wavelengths;  '  according  to  a  recent  theoretical  treatment, 

4 

this  may  very  likely  be  due  to  the  effects  of  the  finite  outer  scale  of  turbulence. 

Similar  data  are  given  in  Figure  3  for  4880  X.  The  logarithmic  slope 

5 

beyond  saturation  is  near  the  (-1/6)  value  predicted  by  deWolf.  However,  as 
discussed  below  and  in  following  sections,  corrections  are  required  for  beam 
refraction  and  covariance  scale  effects,  and  the  agreement  is  probably  for¬ 
tuitous. 

Since  the  strength  of  turbulence  depends  significantly  upon  the 

height  above  ground,  a  more  meaningful  presentation  of  the  data  is  obtained 

by  correcting  for  beam  refraction  and  earth  curvature.  A  useful  abscissa  is 

2 

then  the  theoretical  variance  cr^,  as  predicted  from  the  Rytov  or  first-order 
theory.  *  If  we  define  the  path-length  variable  as  (x),  we  may  write ^ 


-4- 

/ 


0-  2  =  0.  56  k7//6 
T 


i 


C  2M(f)5/b 

n  L 


(L-x) 


5/ 6 


dx 


(1) 


where  L  is  the  total  pathlength  and  k  is  the  optical/infrared  wavenumber. 

7  2 

In  accordance  with  Wyngaard  and  Izumi,  the  height  dependence  of  is 
taken  as  the  minus  4/3  exponent.  We  may  then  account  for  beam  refraction 
and  earth  curvature  by  writing 


2  „  e,  ,  7/6  4/3  _  2,  v  T  -5/6 (■ 

or_  =  0.  56  k  z  C  (z  )  L  * 

T  o  n  o 


yo 


z(x)"4//3  x5/6  (L-x)5//6  dx. 


(2) 


where  z  is  the  beam  height  and  zq  is  a  reference  height. 

It  still  remains  to  determine  z(x),  which  depends  upon  the  vertical 
gradient  in  the  refractive  index,  and  hence  on  the  temperature  gradient.  Our 
measurements  of  the  vertical  temperature  distribution  show  that  the  tempera¬ 
ture  can  be  reasonably  represented  as  linear  over  the  range  involved.  The 
gradient  (B)  is  small,  and  the  resultant  refractive  index  may  also  be  approxi¬ 
mated  as  being  linear  with  height: 


^  B  .  (3) 

dz 


The  trajectory  of  the  beam  through  such  a  medium  is  parabolic.  With  the  use 
of  the  transmitter  and  receiver  heights  as  boundary  conditions,  and  the  addition 
of  the  earth's  curvature,  the  path  may  be  found  for  each  measured  value  of  B. 
Representative  profiles  are  given  in  Table  III.  The  lowest  beam  heights  under 
negative -gradient  conditions,  which  occur  during  the  high  turbulence  part  of  the 
day,  were  comfortably  above  ground  cover  and  agreed  with  visual  observations. 
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TABLE  III.  Beam  Height  (in  Meters)  Vs. 


Temperature  Gradient  and  Distance. 


Thermal 

Gradient 

Distance  from  Transmitter  (km) 

(°K/m) 

0 

1 

2 

3 

4  5 

6 

+  0.4 

6.  1 

7.3 

7.  8 

7.  7 

6.  9  5.  5 

3.5 

-0.  06 

6.  1 

5.  0 

4.  2 

3.  6 

3.3  3.3 

3.  5 

-0.  2 

6. 1 

4.  5 

3.  3 

2.7 

2.5  2.3 

\  5 

ith  the  insertion  of  the  proper  trajectories  into  Eq.  (2),  the  plots 
of  Figures  4  and  5  are  obtained  for  10.  6  um  and  4880  &  respectively.  A 
linear  regression  analysis  for  Figure  5  yields  a  log-log  slope  of  -0.48.  with 
a  correlation  coefficient  of  0.  78,  and  no  apparent  asymptote.  This  does  not 

support  a  recent  theoretical  prediction5  of  (-1/6),  and  there  does  not  appear 
to  be  an  asymptote. 

We  question  the  results  of  Ref.  5  on  the  grounds  that  they  strongly 
involve  the  inner  scale)  in  particular  they  predict  saturation  at  vanishingly 
small  turbulence  levels  for  inner  scales  approaching  zero,  and  this  does  not 
seem  reasonable.  However,  as  pointed  out  in  the  following  sections,  there  is 
an  apparent  contribution  to  the  log  amplitude  variance  at  4880  &  which  is  due 
to  very  small  scintillation  correlation-scales,  and  both  electronic  and  (6  mm) 
receiver  aperture  filtering  have  reduced  this  component  in  the  present  data. 

A  correction  for  this  effect  will  decrease  the  slope  magnitude  in  Figure  5  to 
an  unknown  but  possibly  significant  degree. 

A  plot  of  the  experimental  variance  at  4880  R  vs.  tha*  at  10.  6  is 
shown  in  Figure  6.  Similar  plots  for  shorter  wavelengths  appear  in  Refs.  2  ‘ 

and  8,  and  illustrate  the  failure  of  the  k?/6  dependence  of  Eq.  (1)  in  satura¬ 
tion. 
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B.  Log  Amplitude  Covariance 


the  transverse  log  amplitude  correlation  length  (r  ).  which  is 
defined  from  die  1/e  point  of  the  measured  covariance  curved  is  shown  as  a 
unction  of  Cn  in  Figure  7.  A  significant  effect  is  observed  at  stronger 
turbuience  levels,  where  ra  is  seen  to  markedly  decrease  an  '  increase  at 

.  pm  and  4880  X  wavelengths  respectively.  This  agrees  w  a  trend 
observed  in  earlier  work.  C 

In  order  to  better  understand  the  significance  of  these  results,  it 
IS  necessary  to  consider  the  detailed  shape  of  the  covariance  curves.  Typical 
curves  for  10.  6  pm  are  shown  in  Figure  8.  A.  low  c  *.  the  function  approxi¬ 
mates  the  theoretical  prediction  except  for  large  separations,  where  outer 
scale  effects  may  enter.  However,  as  C,*  increases,  two  new  and  highly 
isparate  scale  lengths  emerge,  such  that  for  strong  turbulence,  the  scintil¬ 
lation  patches  are  ,ui.e  small  with  a  significant  residual  correlation  over 
arge  separations.  This  unusual  manifestation  of  multiple  scattering  was  pre- 
dieted  in  a  recent  two-dimensional  analysis  by  Brow,/  i,  may  be  noted  that  the 
change  in  vertical  beam  trajectory,  such  that  the  highest  turbulence  level  is 
encountered  near  midpath  at  high  overall  turbulence  strengths,  would  imply 

h  emergence  of  a  single  intermediate  scale  sire  and  is  hence  no,  pertinent 
ii©  re* 

This  effect  is  even  more  pronounced  at  4880  1,  as  shown  in  Figure  9 
act  the  effect  is  so  extreme  as  to  suggest  that  significant  spatial  and  elec- 
romc  fUtermg  of  the  small  scintillation  patches  is  occurring  a,  high  turbulence 
levels,  this  is  confirmed  in  the  next  section,  and  indicates  thate.g.  for  curve 
’  a  very  sharp  initial  fall-off  has  been  obscured  and  the  normalisation  is 
incorrect.  This  explains  the  behavior  seen  in  Figure  7.  Visually,  at  high 
turbulence  level,  die  very  large  correlation  patches  are  observed,  and  the 
very  small  patches  are  washed  out  due  to  eye  response.  Note  that  these  results 
,rnp  y  very  poor  receiver-aperture-smoothing,  as  will  be  confirmed  below. 


Scintillation  Spectra 

The  scintillation  power  soectra  at  in  a  ~  i 
and  1 1  far  «,  spectra  at  10.  6  pm  are  shown  in  Figures  10 

e  same  turbulence  levels  as  in  Figure  8  and  ^  k 

x  1“ure  a,  and  corroborate  the 

behavtor  shown  in  ih.  Utter  figure.  As  the  hirbulence  increases.  much 
-re  energy  appears  a,  high  fluencies;  a  concurrent  shift  to  low  frequen- 
ctes  ts  obscured  by  the  lo„-frequency  cutoff  of  the  spectrum  analyzer 

pass  off  ^  rSU'tS  "  4880  *  <FigUre8  12  “*  »>  to  encom- 

enter  F  ^  Wgh  '°  -  —me  scales 

Plot  g‘ng.  lgUre9'  A"  -ception  is  the  transitional  curve  ,B,.  A  similar 
Plot  extending  to  1  KHz  is  not  substantially  different,  and  quantitative  spectral 

results  above  that  fluency  were  not  obtained  due  to  electronic  and  spatial 
(6  mm  receiver)  filtering.  P 

scintilla,-  The  £aC‘  that’  a‘  the“  V"y  ‘arge  “tegrated-path  turbulence  levels, 
scintillations  occur  a,  frequencies  above  one  kilohertz  was  verified  by  direct 

Observation  of  the  signal  (Figure  14,.  This  interesting  behavior  confirms 

effect/  T1  77  CUrVC  <D)  in  FigUr6  9  haS  bee"  °bSCUred  *  tnstrumentai 
to  ts.  A  similar  display  at  10.  6  is  given  in  Figure  15.  The  emerging 

OW  frequency  components  were  also  readily  visible  on  appropriate  time  scales. 

It  is  recognized  that  the  spectral  results  of  Figures  10-13  are  of 

only  semiquantiUtive  signifies  since  the  wind  velocity  was  variable  over 

e  time  periods  spanned.  An  attempt  to  normalize  out  this  effect  as  in  Ref  2 

was  unsuccessful,  due  to  the  complicated  nature  of  the  covariance  and  spell 

curves  and  the  corresponding  inadequacy  of  single-parameter  descriptions 

(e.  g.  r  ). 
a 

D.  Probability  Distributions 

Earlier  theoretical  treatments  predicted  Rqvip;0i,  j 

at  large  turbulence  strenzth.  *  r  •  8  amplitude  statistics 

thatth  •  .  e"Ce  S‘reneths-  E*Psrunental  values  (Figure  16,  have  confirmed 

scintillations  remain  substantially  log  normal,  with  no  indication  of 

2 Tthe  RaT-  hThiS  aSreeS  WUh  m0r'  reCent  'reatmentS-  5  “  “T  be  noted 
that,  in  the  Rayleigh  case,  the  ■'dynamic  range"  between  two  probability  levels 


*  pmr  '• 


is  an  invariant  number  and  that  the  present  data  are  completely  sufficient  to 
rule  out  a  Rayleigh  distribution. 

E.  Receiver  Aperture  Smoothing 


The  covariance  results  of  Figures  8  and  9  suggest  that  receiver 
aperture  smoothing  will  be  very  poor  for  high  integrated -path  turbulence 

levels.  This  was  directly  confirmed  with  large-receiver  measurements,  as 
shown  in  Table  IV. 


TABLE  IV.  Typical  Receiver  Aperture  Smoothing  Results  at  4880  R 

The  receiver  smoothing  factor  R  is  the  log  amplitude  variance  for 
a  32  cm  receiver,  normalized  by  that  for  a  small  (6  mm)  receiver. 


0.46 

0.35 


5.4  x  10 


0.70 


III.  Transmitter  Aperture  (Finite  Beam)  Effects 
A.  Introduction 


The  effects  of  turbulence  on  finite  beams  include  beam  wander, 
instantaneous  beam  spread,  scintillations,  and  wavefront  distortion.  Wander 
and  spread  affect  mean  target-irradiance,  and  wander  and  scintillations 
cause  fading.  Unfortunately,  these  phenomena  have  not  been  well  understood, 
and  theoretical  predictions  in  the  literature  are  characterized  by  apparent 
inconsistencies  and  discrepancies--^  some  cases  due  to  inadequate  definition 
of  the  conditions  under  which  they  apply. 

In  this  section,  we  will  attempt  to  remove  the  discrepancies  and 
unify  the  predictions.  We  will  then  define  a  related  experimental  program. 

B.  Reciprocity 

It  is  physically  apparent  that  beam  wander  is  a  geometric -optics 
phenomenon,  while  instantaneous  spread  may  be  due  to  either  diffraction  or 
multiple  refraction  by  the  turbulence.  Let  us  consider  a  transmitter  optic 
and  target  point,  as  in  Figure  17a;  the  conceptual,  reciprocal  case,-10  involv¬ 
ing  the  point  as  a  coherent  source  and  the  transmitter  as  an  imaging  or 
optical  heterodyne  receiver  system,  is  shown  in  Figure  17b.  It  is  a  conse¬ 
quence  of  reciprocity  that  beam  wander  and  spread  in  the  target  illumination 

system  are  respectively  related  to  image  dancing  and  spread  in  the  reciprocal 
heterodyne. 

In  particular,  if  the  beam  is  nearly  diffraction  limited  but  wanders 
off  target  due  to  atmospheric  effects,  a  virtual  target-point  at  the  new  beam- 
center  will  yield  a  diffraction  limited,  centered  image  (or  high  heterodyne- 
efficiency)  in  the  reciprocal  system;  it  is  clearly  implied  that  the  original 
target  point  will  yield  a  diffraction  limited,  non-centered  image.  This  sug¬ 
gests  a  means  of  eliminating  beam  wander  through  tracking,  as  discussed  in 
later  sections. 

The  beam-image  reciprocity  enables  us  to  make  use  of  theoretical 
descriptions  of  image  behavior  in  order  to  investigate  finite-beam  wander  and 
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spread  at  the  target.  The  application  of  reciprocity  will  also  aid  in  the  under- 
Shtndrng  of  target  scintillations  vs.  transmitter  apermre  diameter,  since  the 
P  enomena  reciprocal  to  scintillation  include  "atmospheric  modulation  noise" 
or  coherent  fading  in  an  optical  heterodyne,  and  this  has  been  analyzed. 

c.  Parameters  and  Structure  Functions 

The  basic  physical  parameters  in  this  discussion  are: 

inner  scale  of  turbulence  ^ 

optical  diameter  ,° 

b 

strength  of  turbulence  c  2 

pathlength  (assumed  horizontal)  ^ 

optical  wavenumber  ^ 

Important  quantities  defined  on  these  parameters'  include  the  wave  structure 

r1:;  DIPl'  lh<!  Pha"  Slr“CtUre  'UnCtt°n  V  <=  ’•  -<*  >  coherence  scale  p 
a  1  of  which  apply  at  the  optic  when  (conceptually,  Inching  at  the  point  source° 

O  more  17b.  The  (asymptotic,  behavior  of  finite  beams  will  depend  upon 
he  relative  values  of  and  b.  i.  e.  upon  which  of  the  si*  possible  inequali. 

tes  between  these  parameters  applies  in  any  given  case.  The  corresponding 
structure  functions  and  P[)  values  are  summarized  in  Table  V. 

With  regard  to  these  structure  functions,  i,  should  be  pointed  out 
*at  Lutomirski  and  Yura  have  calculated  corrections  due  to  outer  scale 
'  effects,  which  they  have  shown  to  be  significant  for  p  >  0. 1  L  m  par 
ticular,  a.  this  value  of  p  the  exponent  is  closer  to  3/2  thin  5/3.  “and  since 
stmple  image  dancing  or  beam  wander  corresponds  t  exponent  of  2.  outer 
scale  effects  may  in  some  cases  lessen  the  advantage  of  wander-tracking.  Also 
L  is  no,  appreciably  greater  than  a  critical  distance  which  is  outer-scale 
dependent,  then  the  value  of  p,.  die  atmospheric  MTF  at  large  p.  and  the  he, ere 
dyne  efficiency  (or  image  resolution,  at  large  aperture,  are  all  increased  over 
values  computea  without  this  outer-scale  correction. 
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TAB.J:  V.  Structure  Functions  and  Cohe 


rence  Scales  for  Point  Sources 

Note:  The  coherence  scale  Pq  is  defined  such  that  D(p  )  =  2. 

The  pertinent  value  of  separation  ( p  )  in  D(p)  is  p  =  b,  the  optical 
diameter. 


WAVE  STRUCTURE  FUNCTION  D 

( 1 )  p  >  b  >  f 
o  o 

b  *  po  >(o 

D(b)  =  1. 1  k2  L  C  2  b5,/3  =  2 

P0=  (0.545  k2  LG  2)'3/5  .  !.44  Lc  2  -3/5 

n 

(2)  p  <  b  <  l 
o  o 

b  <  p  <  t 

o  o 

D(b)  =  0.  62  L  k2  C  2  b2  1  "1//3s  2  f~\Z 

“  ° 

PQ  =  (0.  309  Lk2C  2f  -1/3)-l/2  _.  U80(Lk2c  2|  -1/3  - 1/2 

n  o 

(3)  p  <  i  <  b 

o  o 

b  <  i  <  p 
o  ro 


(4) 

(5) 

(6) 

(7) 


Structure^  Son  not  readily  available  in  literature.  For  complete 
ness,  this  will  Be  further  investigated  at  a  later  date. 

PHASE  STRUCTURE  FUNCTION  D 

4> 

D<f>  =  D  when  b  »  (L/k)l/2  (near  field) 

(c 

rj*  =  f  when  b  «  (Wk)1/Z  (far  field) 
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These  relations  ap  >ly  directly  to  case  (1)  above,  and  we  surmise 
that  they  apply  to  case  (2).  This  will  be  further  investigated  at  a  later  date. 

NOTE:  Case  (1)  obviously  has  the  greatest  practical  importance. 


D.  Beam  Wander  and  Spread 

There  are  a  number  of  alternative  definitions  of  beam  wander  and 

instantaneous  spread.  1  If  tv»  divide  by  the  path  length  in  order  to  use  an 

angular  description,  the  spread  definitions  may  alternatively  involve  the 

second  moment,  half-power  width,  or  on-axis  irradiance;  similarly,  the 

wander  definitions  include  the  second  moment,  and  image-dance  or  ar.gle-of- 

17 

arrival  designations  for  the  reciprocal  system.  In  most  cases,  the  alterna¬ 
tive  definitions  yield  similar  results  which  differ  only  by  numerical  constants. 

In  the  present  discussion,  we  will  invoke  the  simplest  possible 
definitions  in  order  to  clarify  breakpoints  and  asymptotic  dependencies;  we 
believe  that  these  results  are  dimensionally  correct,  and  that  any  necessary 
numerical  constants  may  be  readily  added  at  a  later  date. 

We  point  out  that,  regardless  of  the  definitions  used,  the  mean 

square  total  spread  angle  is  the  sum  cf  the  mean  square  wander  and  instantaneous 
spread  angles. 

Wander 


The  simplest  definition  of  image  dancing  or  wavefront  angle-of- 
arrival  (tilt)  variations  is^' 


*2  = 


kV 


(10) 


Invoking  the  reciprocity  principle,  we  also  utilize  this  definition  to  describe 
beam  wander.  It  will  be  seen  that  the  results  are  dimensionally  consistent 
with  other  treatments  in  the  literature,  15’  1  6'  18’  20“23  and  they  will  clarify 
regions  of  validity  of  those  treatments. 
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A  summary  of  the  mean-square  angle  subtended  by  beam  wander, 
as  defined  by  Eq.  (10),  is  given  in  Table  VI.  This  geometrical  optics  effect 

is  wavelength  independent,  and  depending  upon  the  relative  t  izes  of  b  and  |  , 

.  -1/3  -1/3  ° 

is  proportional  to  b  or  1 

o 


TABLE  VI.  Mean  Square  Eeam-Wander  Angle. 


(1)  b  >  i 


i  <  b  <  p 
o  ro 


$2  =  :.  1  L  C  2  b"1^3 
n 


b  »  (L/k)  '  (near  field) 


=  0.  55  L  C  2  b"1/3  b  «  (L/k)1/2  (far  field) 


f  <  p  <  b 
o  ro 


Same  as  above,  but  beam  spread  may  predominate 
(see  text). 


p  <  I  <■  b 
o  o 


No  structure  function  immediately  available;  beam- 
spread  predominates. 


(2)  b<| 


b  <  p  <  f 
o  o 


4>2=  0.  62  L  C  2  f  _1/3 
n  o 


(near  field) 


•  * 


o 


p  <  b  <  £ 
o 

Same  as  above  (13) --multiple  refractive  effects 
(see  beamspread). 

b  <  t  <  p 
o  ro 

No  structure  function  immediately  available; 
geometrical  optics  gives  above  result  (13). 


Instantaneous  Spread 

For  b>p^,  the  transmitter  aperture  encompasses  a  number  of 
(reciprocal)  coherence  diameters,  and  diffraction  from  one  such  diameter 
yields  the  spread: 


The  results  are  summarized  in  Table  VII.  Note  that  for  p  <  b  <  |  the 

o  o’ 

multiple  refraction  from  the  inner  scale  is  indicated,  and  the  functional 
dependence  is  identical  to  that  for  beam  wander. 


TABLE  VII.  Mean  Square  Instantaneous  Beam-Spread  Angle  (b  > 


po) 


b  >  p  >  l 
o  o 


0.48k2/5L6/5C  12/5 

u  n 


(15) 


1  >  b  > 


=  0.  31  L  C 


2(  -1/3 


n  o 


(16) 


b  >  /  >  p 

o  1  o 

No  structure  function  immediately  available. 
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For  b  <  po,  there  is  no  simple  analytical  result  for  instantaneous 
spread.  ^However,  for  f  o  <  b  <  P(),  the  result  is  implied  in  heterodyne  tilt¬ 
tracking  or  short-exposure  imaging25  analyses;  the  reduction  in  heterodyne 
efficiency  or  image  resolution  which  remains  after  instantaneous  image  dancing 
(beam  wander)  is  removed  is  a  direct  indication  of  image  smearing  (beam 
spread),  to  particular,  as  the  aperture  is  increased,  the  higher-order  or 
nonlinear  wavefront-distortion  terms  (i.  e.  in  the  reciprocal  system)  become 
more  important,  and  spreading  increases  while  wander  decreases  in  accord¬ 
ance  with  Eq.  (11).  to  general,  in  this  regime  the  short-term  spread  is 

smaller  than  either  (1/kpJ2  or  the  wander  angle.  This  will  be  discussed 
further  in  a  later  section. 

Total  Spread,  Breakpoints.  Asymptotic  Behavior 

The  mean  square  long  term  spread  is  defined  as 


Q'Z  =  ez  +  *2 


(17) 


If  we  consider  the  most  important  practical  case,  i.  e.  negligible  inner  scale 
effects  (f  o  <  po,  b),  we  may  clarify  the  discussion  by  constructing  an  asymp¬ 
totic  diagram  of  total  spread  vs.  b,  showing  breakpoints  between  the  various 
regimes.  If  we  add  to  the  previous  definitions  the  mean  square  aperture  dif¬ 
fraction  angle  (kb)  ,  and  assume  the  more  interesting  case  of  the  near  field, 
we  may  use  Eqs.  (11)  and  (15)  to  obtain  Figure  15.  For  the  particular  numerical 

constants  used  here  (i.  e.  unity  in  Eqs.  (10)  and  (14)  ),  the  breakpoint,  are  given 
by 


b.  =  0. 65  p 
1  r  o 

b,  =  8.4  p 

Vbl  =  12-’ 


(18). 


.-16- 


> 


»  • 


Hence  there  is  always  a  wander-predominating  regime,  with  an 
tipper-lower  breakpoint  ratio  which  is  independent  of  all  parameters.  Note 
that  for  b  »  po,  the  spread  is  entirely  determined  by  the  atmosphere,  and 
ts  aperture-independent;  this  may  be  termed  "aperture  saturation. " 

Within  this  wander-predominating  regime,  there  exists  the  possi¬ 
bility  of  reducing  long-term  spread  by  cancelling-out  beam  wander.  This 
cancellation  may  be  achieved  by  making  the  outgoing  transmitter  wave  angle 
(Fig.  17a,  correspond  to  the  incoming  angle  from  the  target  (17b).  A  reflec¬ 
tion  off  the  target  may  be  utilized,  and  in  principle  the  targe,  need  not  be 
cooperative.  If  such  a  tracking  system  is  used,  then  beam  wander  is  elimi¬ 
nated  and  the  broken  lines  in  Figure  18  apply.  There  is  then  an  optimum 
aperture  lor  minimum  long-term  spread,  and  we  note  that  an  Inversion  of  the 
diagram  of  Figure  18  is  consistent  with  plots  of  tilt-tracking  heterodyne 
signal  power  or  short-term  image  resolution?5  This  will  be  further  dis¬ 
cussed  below,  in  the  context  of  mean-target  irradianco. 

It  is  of  interest  also  to  consider  the  far-field  case.  If  b  <  (L/k)1^2 
the  phase  structure  function  is  reduced  to  half  the  total  wave  structure  function, 
the  other  half  being  comprised  of  the  amplitude  structure  function,  implying 
substantial  scintillations.  The  beam  wander  is  reduced  by  a  factor  of  two. 
and  the  breakpoints  corresponding  to  (18)  are 

b  =  0.  99  p 

1  o 

b,  =  1.04  p 

o 

b./'b  r  1.05 

£  I  / 1  n\ 


Hence,  the  wander  regime  essentially  disappeared,  and  instantaneous  beam 
spread  is  practically  a  nil  consideration. 

8  7/0  UA  ‘hat  ‘he  reC>ui''ement  P0  »  2  *  equivalent  to  requiring 

n  k  L  <e  1,  i.  e.  that  the  amplitude  scintillations  not  be  saturated. 

We  point  out  that  the  parameter  r^  which  appears  in  analyses  of 
imaging  or  optical  heterodyne  performance  in  turbulence,  is  simply  related  to  p 
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(20) 


2.  15 


*  intersection  of  the  aperture  and  atmospheric  diffraction  asymptotes  in 

Figure  18  occurs  at  b  =  po  for  our  inexact  numerical  constants  and  unspeci- 

>ed  aperture  Urination;  ro  simply  represents  an  exact  value  for  the  particula, 
case  of  a  uniformly-illuminated  aperture. 

The  relationship  (20)  is  true  in  either  the  plane  or  spherical  wave 
cases,  and  we  note  additionally  that 


=  1.8 


Ithough  the  heterodyne  analyses  were  originally  carried  out  for  plane  wave 
sources,  the  results  carry  over  directly  to  the  spherical  wave  (point  source) 
case  which  is  appropriate  for  reciprocity  reasonings;  we  merely  change  the 
va  ueof  po  or  rQ  accordingly.  This  is  also  true  of  the  scintillation  ana.ysis 

:  eCtt:  F  beW-  —  may  he  immediately  extended 

to  a  nonhorizontal  pad,,  by  utilizing  appropriate  structure  functions  and  values 

Of  p  , 

o 

E.  Mean  Target-Irradiance 

The  application  of  reciprocity  tells  us  that  the  average  target-irra- 
>ance  re  one-to-one  related  to  the  image  resolution  or  the  effective  received 
power  P  m  a  heterodyne  (v  SNR, .  In  fact,  we  expect  that  P'1  „  g,\  We 
replot  the  quantity  P,  normalized  by  that  for  an  infinite  aperture,  in  Figure  19 
(a  en  from  Ref.  25),  and  note  the  consistency  with  Figure  18.  In  particular, 

“  tells  us  that  the  improvement  in  near-field  mean  irradiance  which  can  be 
Obtained  with  wander-cancellation  and  an  optimum  aperture  (for  any  given  con- 
itron  or  value  of  P)>,  is  *  6  dB,  relative  to  that  for  a  very  large,  static  aper¬ 
ture.  This  occurs  a,  b  .  3.  8  r<>.  where  the  instantaneous  abnospheric  beam 
spread  equals  the  aperture  diffraction  spread  (Figure  18). 

Note  that  if  po  <  (L/10  1/2.  the  optimum  b  will  not  be  well  in  the 
near  field,  and  the  above  advantage  cannot  be  obtained.  Hence,  equivalently. 


the  amplitude  scintillations  must  be  well  below  saturation  for  wander-cancel- 

lation  to  be  advantageous;  this  may  often  be  of  more  practical  importance  at 

10.  6  pm  than  at  shorter  wavelengths,  and  will  depend  upon  the  other  parameters 
of  the  situation. 

We  note  that  the  atmospheric  spread  and  wander  angles  implicit 
in  the  heterodyne  analyses  and  Figure  19  are  those  deduced  from  mean  on-axis 
irradiance.  This  does  not  exactly  agree  with  our  wander  definition  (Eq.  10) 
Also,  in  the  heterodyne  analyses  the  optimum  tilt-tracking  is  defined  in 'terms 
of  minimizing  the  integral  of  the  ..  -adratic  phase  function  over  the  aperture 
which  is  not  exactly  equal  to  canceling  D  (b).  However,  for  practical  pur- 
poses  the  differences  in  these  definitions  are  small. 


F.  Fading 

In  the  above  discussion,  we  have  considered  angular  wander  and 
spread  and  their  effects  on  mean  target-irradiance.  Fluctuations  in  this  irra¬ 
diance  arise  due  to  both  wander  and  scintillations,  and  are  also  of  great  impor¬ 
tance.  Unfortunately,  these  finite-beam  effects  have  been  poorly  understood, 
and  in  this  section,  we  will  attempt  to  clarify  these  considerations. 


Scintillation 


An  apparent  contradiction  in  the  theoretical  literature  is  the  following 
erturbation  (Rytov)  scintillation  analyses  have  predicted27’29  that,  for  a  focused 
near-field  transmitter  over  a  horiaontal  path,  the  scintillations  decrease  arbi¬ 
trarily  as  the  aperture  is  made  arbitrarily  large.  However,  the  application  of 

rOpinWneaei..  • _  I  111? 


reciprocity  ana  coherence  theory11'  12  results  in  r 

y  results  in  the  prediction  of  an  "atmos 


*  - - 

phene  modulation  noise"  or  cohprpni  fa^inn  v.  l.  • 

2  erent  fading,  which  increases  asymptotically 

as  b  . 

!n  a  recent  report,  30  we  speculated  that  the  true  conditions  for  appli- 


,  - - ^  ‘Vi.  " 

cability  of  the  Rytov  analysis  become  progressively  more  severe  as  the  aperture 

inrroecas.tLi.L.  '  31 


-  ,  - —  w  » w*  v  uic  operu 

increases;  this  has  now  been  directly  verified  by  Gochelashvily.  31  who  shows 
that  such  analyses  are  valid  only  for  D+  (b)  «  1,  i.  e.  b  «  P().  Dimensionally, 
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it  was  shown  in  Ref.  30  that  tins  is  equivalent  to  requiring  that  the  point-source 
or  spherical- wave  log  amplitude  variance  a  due  to  scintillations  satisfy 


0.  124  C  2k 7/6Ln/6  2 


n 


^  <<  constant  X  (transmitter  Fresnel  number) 


Hence,  as  the  turbulence  increases,  the  maximum  aperture  for  which  the  scin¬ 
tillation  reduction  will  be  experienced  will  decrease;  this  is  consistent  with  the 
onset  of  a  further  noise  mechanism  at  larger  apertures,  i.  e.  coherent  fading. 

We  now  express  the  normalized  variance  of  (linear)  irradiance  a  2 


at  the  target  of  Figure  17a  following  Fried's  approximate  heterodyne  noise 

i  •  12 

analysis: 


I 


2  4a  (r 
(r  -  e 
I 


s  f(b/r  ) 
o 


(23) 


where  cr^  "  is  the  point-source  scintillation  above,  a  is  an  ''aperture  smoothing" 
factor  which  represents  the  Rytov  prediction  of  decreased  scintillation  with 
increasing  aperture,  and  f(b/r^)  is  the  coherent  fading  contribution  ( ^b2  for 
b  >:>  rD**  The  situation  is  summarized  in  the  approximate  asymptotic  diagram 
of  Figure  20.  It  is  apparent  that  the  optimum  transmitter  aperture  size  is  of 
the  order  of  r^,  and  that  a  substantial  decrease  in  target  irradiance  fluctuations 
may  be  achieved  providing  r  »  (L/k)1//2.  We  suspect  that  the  indefinite 
increase  in  cr^  which  is  predicted  for  increasing  b  is  incorrect  and  indicates 
a  further  breakdown  of  the  theory  in  the  multiple  scattering  region.  This  is 
intuitively  apparent,  and  is  supported  by  measurements  in  Ref.  30,  where  it 
was  found  that  scintillations  in  this  regime  (b  »  r  )  approached  those  for  a 
point  source. 

It  is  interesting  to  note  that  the  extreme  focus -criticality  which  is 

predicted  and  observed  for  the  smoothing  of  fluctuations  (b  <  r  )  is  con- 

o 

sistentwith  the  reciprocal  heterodyne  viewpoint;  this  criticality  is  equivalent 

to  the  wavefront-matching  requirement  between  the  incoming  and  local  oscillator 
32 

beams.  We  also  notice  that  the  new  (coherent  fading)  noise  mechanism  is  a 
manifestation  of  atmospherically  induced  beam  breakup  (Figure  18). 


(22) 


It  was  pointed  out  earlier  that  the  requirement  r  »  (L/k)^2  is 

G 

equivalent  to  requiring  that  point-source  scintillations  be  well  below  saturation. 
If  that  is  not  the  case,  the  scintillation  reductions  shown  in  Figure  20  cannot  be 
observed;  beam  breakup  will  take  over  before  aperture -smoothing  occurs. 

We  now  argue  that  beam  wander  (wavefront  tilt)  has  been  implicitly 
removed  in  the  analyses  leading  to  Figure  20,  i.  e.  in  Refs.  12,  27-29.  With¬ 
out  this  removal,  a  new  fading  mechanism  due  to  wander  needs  to  be  added; 
the  measurements  in  Ref.  30  indicate  that  the  resultant  logarithmic  fluctuations 
may  be  an  order  of  magnitude  more  severe  than  the  point-source  scintillations  ! 
We  now  consider  this  point  in  more  detail. 

Wander  Fading 

33  34 

35^  3  recen^  analysis>  Titterton  ’  has  extended  an  analysis  by 
Esposito  to  show  that  the  mean  irradiarce  I  and  irradiance  variance  <r  2 
(on  the  long-term  axis)  due  to  wander  are  given  by 

~  _  i 

I  =  1  +  2  A  (24) 

o 


2  4  A2 

0*  =  - 

I  4a  +  1 


where  Iq  Is  at  the  instantaneous  beam  center,  and 


^  _  mean-square  displacement 
radius 


(25) 


(26) 


Note  that  for  large  wander,  I  ^  ^  as  expected. 

If  the  aperture  diffraction  spread  is  greater  than  the  instantaneous 
atmospheric  spread  (Figure  18),  we  have  from  Eq.  (11): 
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A  'v* 


D(b). 


(27) 


1.1LC  “  b 
_  n 


2  .  -1/3 


■  2U2 
k  b 


2  b5^3 


(In  Titterton's  analysis,  the  5/3  exponents 
this  with  Eq.  (25),  we  have 


are  approximated  by  2.  ) 


Combining 


cr 


I 


2 


i-D(b)Z 
4D(b)  +  1 


(28) 


Similarly,  at  large  apertures  02  =  (kp^'2  and  we  haye 
A  ^  2(  p  /b)1^3  =  2.3  (D(b)l'1/5 

(29) 

2  21  O'2/5 

CT  *\j  — 

9.2D'1/5+l  '  (30) 


these  expressions  ,28.  30)  may  be  utilised  in  another  asymptotic 
lagram  (Figure  21).  where  the  peak  value  of  tr/ft)  is  a  fixed  number  independent 

m  .Tb  1  T'  “  Wand<!r  CarU:<!Uati0n  iS  "0t  ^Ployed,  this  fading  mechanism 
be  added  to  those  shown  in  Figure  20.  as  crudely  indicated  there.  The 

relative  height  of  the  wander  fading  wil,  depend  on  ,  measurements  have 

s  own  an  order  of  magnitude  increase  in  the  logarithmic  fluctuations  due  to 
this  mechanism.  In  future  efforts  ,„ui  „  1 

..  ettorts,  we  will  seek  more  precise  quantitative  pre- 

lctions  as  a  function  of  the  independent  parameters. 

Summary 


in  summary,  it  is  apparent  that  the  successful  cancellation  of  beam 
wander  has  two  potential  advantages: 

)  The  mean  target  ir radiance  may  be  improved  up  to  6  dB 
over  that  for  a  very  large  aperture 
(2)  The  fluctuations  in  target  irradiance  may  be  very  sub¬ 
stantially  reduced  relative  to  those  for  a  large  aperture. 


°rdCr  '°  aChicv°  *"•  <•»  transmitter  aperture  must  be  variable 

tn  accordance  with  the  value  oi  P().  so  that  a  knowledge  of  the  atmospheric  MTF 
is  critical. 

H.  Further  Comments  on  the  Literature 

There  have  unfortunately  not  been  many  measurements  of  wander  and 
spread,  anyone  of  scintillations  with  wander-cancellation.  Aicaraz  and 
Livingston  have  reported  a  rough  corroboration  of  the  range-dependence  in 
Eq.  Ill  I.  but  a  poor  correlation  with  C^.  It  is  to  be  expected  that  substantial 
contributions  to  beam  wander  will  occur  due  to  turbulence  scales  larger  than  the 

outer  scale,  where  Cn  is  not  a  meaningful  measure.  This  is  especially  true 
for  longer-term  components. 

Kunger  9  has  described  a  modulation  phase  technique  for  measuring 
wander,  and  presents  some  data.  Hansen  and  Madhu40  have  observed  and 
attempted  to  explain  a  curious  effect  involving  large  image  dancing  for  a  retro 
return  relative  to  a  nearby  incoherent  source!  we  have  not  been  able  to  dupli¬ 
cate  this  effect,  and  we  do  not  accept  the  explanation  given.  This  possible 
anomaly  will  be  investigated  further. 


ducted  at  NRL 


c°mprehensive  wander  and  spread  measurements  have  been  con- 
’  ’  but-  unfortunately,  no  attempt  was  made  to  identify 


the  data  in  terms  of  realms  of  I,  q,  b.  pj,  so  that  interpretation  is  difficult. 
For  instance,  data  a,  6328*  |b  »  p<),  are  intermixed  with  data  a,  10.6  pm 
(b  <  pQl.  A  high  degree  of  correlation  was  found  between  wander  at  the  two 
wavelengths,  which  supports  the  geometric,!  view  of  that  phenomenon.  We 
believe  that  these  very  substantial  data  can  be  more  usefully  interpreted  by 

relating  them  to  p  and  (if  known)  f  . 

o  o 

I.  Physical  Viewpoint  of  Transmitter  Aperture  Smoothing  of  Scintillations 

In  this  section,  we  attempt  to  lend  some  insight  into  the  phenomenon 
Of  transmitter-aperture  smoothing  of  scintillations  lb  <  rQ  in  Figure  20).  using 
a  heuristic  physical  argument  based  on  ray  optics  and  interference.  The  degree 
of  scintillation  a,  any  point  is  determined  by  the  percentage  of  incident  rays 
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which  have  relative  phase  shifts  >  tt. 


(1)  Within  the  Focal  Spot-Horizontal  Case 

Consider  a  near-field  source  S  focused  on  a  receiver  plane 
R  at  distance  L,  as  in  Figure  22.  We  know  that  the  physical  optics  solution 
in  a  vacuum  leads  to  a  finite  focal  spot  at  R,  and  we  consider  all  of  the  energy, 
including  this  diffraction  spread,  as  represented  by  rays.  One  such  (arbitrary) 
ray  is  shown  in  the  figure,  and  we  note  that  the  size  of  R  is  determined  by 
the  requirement  that  rays  incident  at  an  angle  equal  to  the  subtense  of  S  have 
substantially  equal  lengths  or  phases  for  all  points  on  R  (in  particular,  at 
both  edges  of  R). 

We  model  the  turbulence  as  consisting  of  small,  weak  refractive 
scatterers  giving  rise  to  scattered  rays  such  as  the  one  shown,  where  the 
(small)  scattering  angle  is  (3,  and  z'  is  the  distance  from  the  scatterer  to  R. 

In  order  to  ensure  that  the  scattered  ray  cannot  interfere  with  the 
nonscattered  rays,  we  require 


k  p  z'  «  1 


(31 


Since  (0  z')  =  R 

max 


kS 


,  we  rewrite  this  as 


kPR  =  «  1.  (5«| 


(32) 


This  simply  manifests  the  way  in  which  R  was  established;  all  non-interfering 
scattered  rays  subtend  S.  Also,  we  note  that  interference  can  occur  only  for 

z'/L  <  I../kS2  (33) 

'V  ' 

Hence,  as  the  source  Fresnel  number  is  increased,  rays  which  can 
interfere  are  restricted  to  larger  scattering  angles  and  a  smaller  scattering 
region.  Thus  for  large  source  apertures,  scintillations  in  the  focused  spot  are 
drastically  reduced. 


(2)  Within  the  (Far  Field)  Spot- -Vertical  Case 

Consider  a  ray  which  is  scattered  near  a  collimated  trans¬ 
mitter,  as  shown  in  Figure  23.  The  scattering  is  assumed  to  occur  through 
an  angle  p  at  a  distance  z"  from  the  source;  since  the  far-fieid  is  by  -’efinition 
angularly- re  solved,  we  compare  the  scattered  ray  with  an  unscattered  ray  to 
which  it  is  parallel.  To  ensure  no  interference,  we  require 

k  ft2  z"  «  1  •  (34) 

SL»ce  Pmax  -  R/  L  within  the  far  field  spot,  we  rewrite  this  as 


k  z" 


k2s2 


«  1. 


«  1 


(35) 


Hence,  scattering  which  occurs  within  the  near  field  of  S  will  cause  no  appre- 

*1  f\ 

ciable  scintillation,  in  agreement  with  theory. 

(3)  Outside  the  Focal  Spot,  at  Focal  Plane 

In  this  region  all  rays  are  the  result  of  scattering,  and  we 
expect  a  high  degree  of  randomization  and  interference,  which  is  observed.  30 
Since  the  region  within  R  scintillates  less  as  S  is  increased,  by  conservation 
of  energy  there  is  then  less  irradiance  appearing  outside  the  spot. 

(4)  Outside  Near  Field  of  Focal  Region 

Consider  a  ray  which  is  scattered  by  a  small  angle  (S/L»  p» 
R/L)  and  bound  for  a  point  P  outside  the  focal  spot,  with  a  correspondingly 
significant  phase  shift  and  ability  to  interfere  (Figure  24).  Outside  of  the  near 
field  of  the  focal  spot,  it  can  overlap  the  unperturbed  rays.  Hence,  for  small 
transmitter  defocusing,  scintillations  on  the  target  become  very  substantial.2^'30 

(5)  Multiple  Scattering 

The  above  discussion  assumes  single  scattering  events,  in  which 
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each  ray  is  scattered  only  once  and  the  basic  field  is  substantially  unper¬ 
turbed.  If  the  scattering  is  stronger,  much  of  the  energy  is  taken  out  of  R, 
and  the  spot  is  spread  with  considerable  opportunity  for  interference  every¬ 
where. 

In  this  situation,  it  is  not  clear  that  the  scintillating  patch  size  is 
the  usual  L/kS,  although  measurements  support  this  covariance  scale  under 
these  conditions.  '  For  an  atmospherically  spread  but  substantially  confined 
beam,  conservation  of  energy  requires  that  the  covariance  must  integiate  to 
zero  over  the  cross-section.  Thus  the  patch  size  cannot  increase  indefinitely 
with  increasing  turbulence  as  it  does  in  Section  II  for  a  point  source.  When 
the  (kp  )’2  atmospheric  spread  predominates,  a  larger  transmitter  may  be 
expected  to  result  in  added  independence  and  decreased  patch  size  (perhaps 


Note  that  as  S  is  increased  and  unperturbed  R  is  decreased,  this 
mechanism  of  substantial  scattering  outside  of  R  is  effective  at  smaller 
scattering  angles.  A  measure  of  the  effect  is  the  value  of  D^(S),  as  per  previous 
discussions  of  the  breakdown  of  Rytov  analyses. 

(6)  Other  Remarks 

Beam  wander  can  be  represented  in  this  model  as  occurring 
due  to  larger  refractive  wedges.  This  effect  does  not  alter  the  preceding  con¬ 
siderations. 

Aberrations  in  the  optical  system  produce  unperturbed  rays  out  of 
R  (Figure  25),  which  can  be  scattered  back  into  R  to  interfere  and  cause 
scintillations,  even  in  the  single-scatter  regime.  It  may  be  shown  that,  for 
scattering  e.  g.  near  midpath,  interference  will  occur  for  scatter**  ancles 

P  >  RS 
L 

An  alternative  and  more  rigorous  approach  to  these  heuristic  argu¬ 
ments  would  be  to  apply  a  Fourier  optical  approach,  whereby  the  turbulence 
is  represented  in  terms  of  thin  phase  screens  having  single  spatial  frequen¬ 
cies,  36  each  of  which  scatters  part  of  each  incident  plane  wave  component 
into  two  new  (weak)  angular  components.  Alternatively,  each  can  be  considered 


to  scatter  part  of  the  unperturbed  gaussian-mode  energy  through  these  two 

angles,  and  the  above  arguments  can  be  applied  in  order  to  determine  the  pos- 
sibility  for  interference  or  scintillation.  ^ 

J*  Experimental  Program 

We  are  now  instrumented  to  demonstrate  the  cancellation  of  beam 
wander  through  the  use  of  a  tracker  with  a  small  targe,  retroreflector.  Using 
the  reciprocity  philosophy,  the  instantaneous,  mean  angle  of  arrival  of  the 
reflected  energy  is  matched  by  the  outgoing  transmitter  wavefront.  We  have 

already  qualitatively  observed  that  this  approach  maintains  the  beam  centroid 
on  the  target  retro. 

In  order  to  quantify  the  results  and  to  explore  the  relationships 
discussed  in  the  preceding  sections,  we  will  perform  measurements  as  out¬ 
lined  in  Table  VUL  These  experiments  will  be  conducted  a,  6328X.  with  and 
without  wander-canceHa, ion,  with  a  near-field  focused  transmitter. 

TABLE  VIII.  Finite-Beam  Measurements- -Target  Illumination 


Parameters: 


Strength  of  turbulence  (C  2,  p  ) 

n  Ko 

Inner  scale  (f  ) 
o 

Aperture  size  (b) 

Range  (L) 


Measurements:  Microthermal  C  2 

n 

Turbulence  spectrum  for  / 

o 

Target  mean  irradiance 
Target  irradiance  fluctuation  statistics, 
variance,  and  power  spectrum 


K.  Angular  Beam  Dithering 

During  initial  experiments  with  the  tracking  system,  a  loop  insta- 
b.my  resulted  in  a  small,  rapid  (1  kHz)  angular  dithering  of  the  beam.  This 
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resulted  in  a  substantial  washing -out  of  visually  observed  scintillations.  The 

explanation  for  this  is  that  turbulence  structure  near  the  transmitter  results 

in  a  scintillation  pattern  which  dithers  with  the  beam,  and  this  component  of 

the  total  pattern  is  hence  being  dithered  faster  than  eye-response.  The  small 

43 

dither  angle  can  be  achieved  electro-optically  at  e.  g.  microwave  rates,  and 
it  is  suggested  that  this  effect  might  be  useful  since  it  would  translate  much 
of  the  scintillation  power  spectrum  to  GHz  frequencies  and  hence  out  of  the 
necessary  response  range  of  many  atmospheric  systems.  We  will  pursue  this 
further. 


IV.  ihe  Effects  of  Turbulence  Intcrmittency 

As  discussed  in  the  preceding  renort  3  Mn0 

°r  n“rm — - — wHicb'u  often  z::::i;r"ncy 

Z  ‘°  lmVe  SUbStan,ia'  a"d  hithe"°  -  -ec,  on  -scintillations 

-  narrow  re6io„s  a,o„s  the 

variations  in  scintillation  ievels  due  to  th  ^  Sta,iS‘iCa‘ 

instantaneous  pad, -realisation  and  ^  re‘ali°"SMP  betW- 

The  usual  theoretical  exnr  •  0"8'term  °r  eDS,!mble  average  conditions, 
in  terms  of  iong .,cr ‘ 7 ’  ^  «<*  turbulence 

certainly  mash  short-terTdTep  777* ^  ^  inValidat,!d  *  ^ 

—  -  -read  ebich  is  ^ 

that  it  '  ‘  diSCU8Si0n>  ^  ““  «».  elite,  and  show 

1  s  ;  ;ra:y  not  be  ~  - 

win  involve  ^  ^  ^  ^ther  worle 

re, teal  and  expertmental  investigation  of  this  latter  topic. 

A-  Propagation  Through  a  Turbulent  Slab  -  Inner  Scale  Effects 

-bs.  coi::“;:rre  in,ermittency  u  ^  * — 

oathlength  L  ■  y  p  opagating  a  quasi-plane -wave  over  a 

from  the  receiver.  The  turbulence  within  the  .1  k 
is  statistically  uniform  with  a  streneth  c  2  „ 

model:1  n  ’  fltS  the  KoI™g°rov  spectral 


-  0.  033  C  V11/3  -*2/k 

n  n  e  rn 


(36) 


where  ^  is  the  refractive  index 
km  ~  5‘  92//*  o*  Ihe  loB  amplitud 


,r  =  2  TT 


‘H 


spectrum,  K  is  the  spatial  wavenumber,  and 
e  variance  of  the  received  beam  is  given  by1'44 


•F  (  K )  4>  (  K  )  K  d 
n 


k. 


(37) 
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*  ^  f+Bi~  t  l«-  Mtl 


WhCre  44^  18  ^  aPPr°Priate  filt«r  function  for  the  path  under  consideration. 
DeWolf  has  given  F(K)  for  the  path  of  Figure  26  as 


F(k)  =  1  -  sin  d)  *2/k]  -  sin  [bK2/k] 

K^d/k 


(38) 


Using  (38)  and  (36)  with  (37),  we  have 


2 

<r 


=  2tt2(0.033)C  2k2d 
n 


r 

dldK  K 

Jo 


-8/3e-K2fo2/(5.  92)* 


sin  f(b  +  d)  K^/k]-8in[bK^/l< 
«2d/k 


(39) 


In  order  to  check  for  consistency  with  well  known  results,  we  let 
b  —  O  and  d  -  L  in  Eq.  (39): 


2 

<r 


(40) 


This  is  equivalent  to  Eq.  (Tl)  in  Ref.  1. 

In  order  to  carry  out  the  integration  in  Eq.  (39),  we  let  y  =  K2 : 


(41) 


This  yields 
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/ 


f  •'*-#••**■* 


Noting  that  f~(~  =  -yy-  f  (-5/6),  this  becomes 


(43) 
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2 


or  b  «  '•and  b/k  «  (  ^  this  expression  reduces  to  (T4)  and  (T5)  of  Ref.  1. 
for  d/kf  ^  | <<J  1  respectively. 

We  are  most  interested  in  a  fairly  small  value  of  d  at  various 
values  of  b.  and  the  extent  to  which  the  inner  scale  affects  results.  We 
arbitrarily  chose  k  =  1.2,  loVV  .  5200  X).  and  let,  range  from  3  to 
30  mm,  such  that  the  corresponding  values  of  ,  \  range  from  108  m  to  10.  8 
km.  The  results  are  given  in  Table  IX.  We  note  that,  for  d  -  11m,  the  inner 
scale  effects  are  substantial  at  b  =  100m  but  much  less  at  b  =  1  km. 

Since  the  scintillations  are  strongest  for  slabs  furthest  from  the 
receiver,  where  inner  scale  effects  are  minimal,  it  may  be  expected  that 
the  average  effects  of  the  inner  scale  over  a  long  pad,  are  not  great.  W.  now 

derm  a  good  approximation  to  Eq.  (43)  which  is  valid  for  most  cases  of  inter- 
est  anc  for  which  the  inner  scale  drops  out. 


TABLE  IX.  Values  of  cr2/^2  from  Eq.  (43).  k  is  taken 


as  1,2  x  10^  m  \ 


b(m) 

d(m) 

i  0(cm) 

<r2/C  2  (m2/3 

100 

11 

0.3 

n 

4. 9  x  1010 

0.9 

3.  1 

1.2 

2.  3 

2.  1 

9.  0  x  109 

5C0 

11 

3.  0 

0.3 

3*4  11 

2.  0  x  10 

0.9 

1.  7 

1.2 

1.  5 

2.  1 

1.0 

1000 

11 

3.0 

0.  3 

7.0  x  loj° 

3.  6  x  10 

0.9 

3.  3 

1.2 

3.  0 

2.  1 

2.4 

1000 

100 

3.  0 

0.  3 

J*9  12 
3.4  x  10i2 

0.9 

3.  1 

1.2 

2.  9 

2.  1 

2.  3 

3.0 

1.  9 
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Referring  to  Eq.  (43),  let  us  assume  that  b  »i  q  2k.  We  thus  have 

5/6  r  11/6 


^  =  0.326  c/kVf-fJJp 


6  k 


r. 


5.  92 


11  d 


b_+_d  j 
k 


•  /  Uw\ 

l  »/ 


L  J 


11/6 


0.31  Cn2k?/6  (b+d)11/fc  -bll/6j. 


D 


As  a  further  approximation,  we  may  require  b  »  d,  and 


write 


o-2  =  0.  57  C  2  d  k7^6  h5/6 


n 


(4 


Which  is  the  central  result  o£  this  section.  We  note  that  the  dependence  on  the 
inner  scale  has  vanished  in  this  approximation,  even  though  d  may  be  <  ,  2k 
and  the  dependence  on  b  is  clear,  fa  Table  X,  „e  compare  results  calculated 

from  Eq.  (45)  with  those  from  the  complete  expression  (43).  The  approximatio 
is  seen  to  be  good  even  for  b  ~  £  2k. 


W 


may  WrUe  E<>-  <45)  in  a  more  form  by  recognising  that  the 


average  Cn  over  the  path  (equal  to  the  long-term  average  measured  at  a  point) 


IS 


C  *  d/L 
n  ' 


(46 


where  L  is  the  total  path  length.  We  thus  have 
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,2=  0. 

57  C  2  k?//6  b5/6  L 

(47) 

n  • 

TABLE  X. 

• 

2  2 

Values  of  <r  /Cn  from  Eq 

s.  (45)  and  (43). 

k  is  taken  as  1.2x10* 

and  i  q  as  3  mm,  such  that  k£  2  =  108m. 

b(m) 

d(m)  c 

r2/c/  to2/3, 

<r2/Cn2  (m2/3) 

(Eq.  43) 

(Eq.  45) 

100 

1 

11 

51 

101 

200 

4.  3  x  109 

4-9x10 

2.7  X  10  ; 
6.3x1°" 

1.6  x  10 

4.  8  x  109 
5.3x10  0 
2.5x10 
4.9x10 

9.  6  x  10 

500 

1 

11 

51 

101 

200 

1.8  x  lo}° 

2.0  x  loJJ 
9.5x10“ 
1.9x10 

4.2  x  10  ^ 

1. 8  x  I0l° 
2.0x10 

9.4x  10  ' 

i.9x  io, ; 

3.  7  x  10  ^ 

1000 

1 

11 

51 

101 

200 

3.  2  x  10  j? 

3.6  x  10 

U7xl°  2 

3.4  x  10  2 

7.0  x  10  ^ 

3.  3  x  l°j? 

3.  6  x  10 

1.7x  10  2 

3. 3  x  1 °i2 

6.  6  x  10 

In  order  to  obtain  the  long-term  average  of  scintillations,  we  calculate  the 

average  of  Eq.  (47)  over  all  values  of  b  such  that  0  <  b  <  L.  This  corresponds 

to  an  equal  probability  for  all  possible  positions  of  the  slab  along  the  path.  The 
result  is 


=  0.  31  C 


7/6  11/6 


(48) 
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w  ich  -  the  well  known  expression  for  a  plane  wave  In  uniform  turbulence  1 
ince  Bq.  (47,  does  not  involve  the  inner  scale,  we  expect  this  result. 

term  a  ^  reSUl'’  ““  '"n°r  h"  °n‘y  3  SmaU  effect  »■>  long- 

.eaa  *»"  *•  turbulence  regions  are  much 

stantial  nfl  °  ’  may  ^  "  f°‘IOW8:  1116  inner  SCale  has  a 

influence  only  for  slabs  near  the  receiver,  but  those  slabs  in  turn 

give  rise  to  a  relatively  small  degree  of  scintillation  and  hence  are  weakly 

weighted  in  the  averaging  process.  The  approximation  is  within  5%  for 

Eo  S  V  =  3m”’  and  d  =  ‘m’  38  determinad  f  >  comparison  of 
q*  Wlth  a  true  Path -averaging  of  Eq.  (43). 

We  thus  conclude  that  the  results  of  long-term  measurements  of 

s  C  will  not  be  appreciably  affected  by  such  intermittent,  although 

uch  a  statement  was  not  obvious  before  detailed  analysis.  However,  there 

7  “  meChanUm  Si8nifiCaat  fluctuations  in  these  qua„- 

,  and  this  will  be  discussed  further  in  the  next  section. 

B.  Further  Definition  of  the  Intermittency  Statistical  Problem 

tr  ^  tte  PreCedi"6  reP°rt  °n  111:8  Pr°gram’  3  We  haya  discussed  the 

‘  °  ^lnlte " tlme “average  measurements  of  scintillations  (log  ampli- 

tu  e  and  microthermai  fluctuatlons  at  a  point  as  ra;dom  p 

3.  We  denote  these  discrete  random  variables  or  measurement  samples 

T  .  Y  u.  order  to  state  the  problem  in  a  completely  abstract  way.  The 
physical  definitions  are 


xt  “  Cn  (T  )  =  const,  x  - 

n  T 


(AT(t)  )2  dt 


yT  -  °"  (T  )  =  const,  x  — 

T 


O 


i2(t)  dt 


(50) 


22 ati  iS  the;vera8tog  ,tee’  *  T<*>  ‘He  differentia.  temperature 

actuations,  and  i  („  denotes  die  log  atnp.itude  scintiUations.  For  r  »  de. 

corre  ation  tunes,  x  and  y  become  normally  distributed. 

*  T:::  7  "*  * ~ — - 

<■)  the  convergence  of  either  one  separately  to  a  mean  can  be 
analyzed  vs.  the  averaging  time  r  and  the  power 
spectrum  or  autocorrelation.  46 

(2)  A  linear  or  nonlinear  regression  analysis  may  be  applied 

betaken  x  and  y,  with  resultant  confidence  intervals 
and  correlation  coefficients.  16,46 

aTfllwT8"'  lhe0rellCal  aDd  eXPerimenlal  eXlenai°n’  -  considerations. 

,3)  What  can  be  said  apriori  concerning  the  statistical  relation- 
ip  between  and  yT  ,  including  data  spread  or  con¬ 
fidence  intervals,  and  correlation  coefficients,  given  a 
theoretical  relationship  between  related  quantities  , 

(4)  How  does  this  relationship  behave  as  a  function  of  averag- 
ing  time  T  ? 

;:an:;::ry’  ^  -  *•  H.evan,  Physical  variables 

needed  for  .meaningful  answers  to  the  above  questions  ? 

above  T°  de'Cribe  ‘he  ‘he0rettCal  relationship  referred  to  in  item  ,3) 

.  ’  aS80Clale  Wlth  ^  Pb“t  quantity  A  T(t)  a  spatial  quantity  A  T(a) 

T  *  18  th'  °Ta'  Pathle"gth  Variabla  n<z<L>.  If  we  assume  the 
ay  or  ypo  esis,  and  the  wind  happens  to  he  along  the  path,  these  two  quan- 

titles  are  simply-related:  * 

AT(t)  =  AT(z/v) 
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where  v  is  the  wind  velocity.  In  general,  however,  the  macroscale  inter - 
mittency  may  not  be  assumed  to  be  isotropic;  that  is,  the  turbulence  regions 
will  not  necessarily  have  equal  dimensions  along  and  perpendicular  to  the 
wind  direction. 

We  now  suppose  that  some  particular  realization  of  the  path  is 
frozen,  i.  e.  ,  that  the  nonuniform  turbulence  field  is  invariant.  From  the 
results  of  the  preceding  section,  we  drop  inner  scale  effects  and  write  for 
a  nonuniform  path:* 


where  f(z)  is  a  weighting  function.  Since  it  is  not  practical  to  know  the  mean- 
square  temperature  fluctuations  along  the  path,  these  expressions  are  not 
directly  useful. 

We  therefore  let  the  overall  turbulence  structure  evolve,  and  assume 
stationarity  in  the  long-term  statistics.  To  the  extent  that  this  evolution  is  slow 
compared  to  the  actual  temperature  and  log  amplitude  fluctuations,  we  may 
utilize  a  finite  t  in  Eqs.  (52,  53),  where  t  is  intermediate  to  the  fluctuation 
and  evolution  time  scales  respectively.  It  would  appear  that,  given  appropriate 
statistics  of  Cn  (r  ,  z)--e.g. ,  the  z -correlation  function  and  temperature  power 
spectrum— the  above  questions  can  be  answered.  The  integral  relationship 
clearly  indicates  less  spread  in  y  than  in  x. 
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C.  Disparate  Time  Scales  in  Random  Processes 


As  a  further  attempt  to  define  the  problems  posed  by  intermittencv, 
we  consider  the  following.  Suppose  that  the  random  process  x(t)  is  made  up 
of  the  product  of  two  processes: 

x(t)  =  x^(t)  x  ^  (t)  ,  (54) 

where  the  frequency  scales  involved  in  the  two  sub-processes  are  highly  dis¬ 
parate,  and  the  higher-frequency  process  has  zero  mean  (Figure  27).  For 
example,  suppose  that  the  power  spectrum  of  x^  ranges  from  10  to  1000  Hz, 
while  that  of  x ^  ranges  from  0  to  10  1  Hz.  Since  the  two  are  multiplicative, 
the  spectrum  of  x(t)  is  at  the  higher  frequencies,  with  low-frequency  side¬ 
bands  on  the  Xj  components  (Figure  28). 

The  power  spectrum  and  autocorrelation  function  of  x  do  not  appre¬ 
ciably  evidence  the  presence  of  the  low-frequency  multiplier,  x^  However, 
averaging -time  considerations  such  as  those  reviewed  in  Ref.  45  will  be 
invalid,  and  the  overall  random  process  appears  nonstationary  in  the  short  term. 
This  breakdown  of  the  methods  of  Ref.  45  is  physically  obvious,  since  any  short¬ 
term  attempt  to  estimate  e.g.  the  mean  of  x2  will  be  invalidated  by  long-term 
trends;  the  formal  breakdown  of  the  theory  is  not  as  immediately  apparent. 

In  order  to  avoid  these  difficulties,  we  may  deal  instead  with  the 
random  variable  x2(t).  The  low  frequency  components  in  x2  will  then  appear 
directly  in  the  power  spectrum  and  autocorrelation  function  (Figures  29  and  30). 
For  example,  for  any  finite-time  estimate  of  x2,  we  may  write 


where  the  additive  terms  ^  and  represent  the  low  and  high  frequency  con¬ 
tributions  to  che  error,  respectively.  These  are  of  course  random  quantities, 
with  f  j  »  f  . 
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The  low-frequency  multiplying  function  x^t)  represents  the  effects 

of  interim  ttency.  In  one  model,  it  consists  of  a  simple  two-level  function.  47 

In  general,  its  probability  distribution  is  of  interest;  the  distribution  of  x 

48  2 

or  AT  is  understood.  This  separation  of  the  temperature  function  into  two 

frequency  or  time  scales  may  be  usefully  applied  to  the  considerations  of  the 

preceding  section.  We  will  of  course  continue  to  refine  our  theoretical  state¬ 
ment  of  this  problem. 

D.  Experimental  Approach 

Our  field  site  has  been  Shortened  to  a  uniform  path  or.  the  order  of 
one  mile  in  length,  and  restored  to  a  point-transmitter,  point-receiver  con¬ 
figuration.  The  microthermal  and  optical  (4880&  and  10.6  pm)  instrumentation 
is  being  adapted  for  digital  computer  processing  of  analog  fm  tape  recordings 
of  the  raw  data  (AT(t)  and  f  (t)  ).  This  facility  is  independent  of  that  being 
utilized  in  the  finite-beam  and  wander-cancelletion  work. 

The  first  experiments  will  consist  of  (x,  y)  regression  determina¬ 
tions  vs.  r,  which  will  be  related  to  the  power  spectra  of  AT  and  |.  Fur¬ 
ther  efforts  may  then  involve  the  measurement  of  the  z-correlation  of  C  2, 
using  two  or  more  microthermal  systems  with  variable  spacings  between  the 
systems,  and  a  second  such  system  is  currently  under  construction.  This 
correlation  measurement  is  equivalent  to  determining  an  experimental  model 
of  the  macroscale  intermittency  of  the  turbulence. 

In  accordance  with  the  preceding  section,  it  may  be  of  interest  to 
treat  the  square  of  AT(t),  to  single  out  the  low-frequency  portion  AT(t)  and 
measure  its  probability  distribution,  and  to  measure  the  statistical  and  spe  ctral 

properties  of  the  squared  quantity.  This  may  be  readily  achieved  on  the  digital 
processor. 
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V.  Publications 


The  following  paper  is  to  be  published: 
on 

- ..  ,m  s— 

Large  Integra  ted -Path  Turbu/ence  ^Awritte^41^^0"  ^e*surements  for 
prepared  for  publication  in  the  near  future.  CrS1°n  S  PapCr  ^  be 
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VII.  Figures 

1.  Experimental  log  amplitude  variance  vs.  time  of  day. 

(X)  488oX  on  .September  14,  1972 

(o)  10.  6|im  on  September  14,  1972 

(.)  488oX  on  August  29,  1972 

The  weather  was  clear  except  for  a  trace  of  ground  fog  from  0600-0700 
on  September  14.  During  the  neutral  period  of  typically  0700-0800,  the 
turbulence  was  poorly  developed. 

2.  Dog  amplitude  variance  at  10.  6  pm  vs.  strength  of  turbulence  at  1.8  m 
height.  The  data  points  include  several  clear  days  of  operation. 

3.  Same  as  Figure  2,  for  4880  X. 

4.  Experimental  vs.  theoretical  (Rytov)  log  amplitude  variance  for  10.6  urn. 

The  abscissa  is  corrected  for  beam  refraction  and  earth  curvature  effects. 

The  line  indicates  the  (cr  2  =  or  2)  condition. 

h,  T 

5.  Same  as  Figure  4,  for  4880  X.  The  equation  of  the  linear  regression  line 
is 


1Og10°E2  B  <-°-22)  -  (0.48)  log1()<rT2 


with  a  correlation  coefficient  of  0.  78. 

6.  Experimental  variances  at  4880  X  vs.  those  at  10.  6  urn.  The  line  repre- 

7 /6  r 

sents  a  k  dependence. 

7.  Transverse  log  amplitude  covariance  length  vs.  strength  of  turbulence. 

The  experimental  and  theoretical  (Rytov)  values  are  for  4880  X  (X, - ) 

and  10.  6  pm  (0, - ). 

8.  Normalized  covariance  curves  for  10.6  pm.  The  normalizing  quantity 
C|  (0)  is  identical  to  cr^,  ,  and  the  l/e  points  are  indicated.  The  broken 

line  represents  the  theoretical  (Rytov)  function. 
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Curve 

c  2 

n 

A 

-15 

4.2  x  10 

B 

3.  6  x  10‘13 

C 

5. 9  x  10‘13 

D 

3.2  x  10"12 

9.  Same  as  Figure  8,  for  4880  &. 


Curve 

A 

B 

C 

D 


n 

4.  5  x  10"15 
3.8  x  10‘13 
6.  3  x  10-13 
3.4  x  10‘12 


10. 


RMS  scintillation  spectra  at  10.  6  jim  for  the  C  2  values 
The  1/e  frequencies  are  indicated.  n 


of  Figure  8. 


11. 


Same  as  Figure  10, 
quency. 


with  the  ordinate  squared  and  weighted  by  the  fre- 


12.  RMS  scintillation  spectra  at  4880  8  for  the  C  2  values  of  Figure  9. 

13.  Same  as  Figure  12,  with  the  ordinate  squared  and  weighted  by  the  fre- 


14. 


Linear  (left)  and  log  (right)  irradiance  scintillations  at  4880  8.  The 
pictures  from  top  to  bottom  represent  C  2  values  A-D  respectively  of 
Figure  9.  The  abscissa  is  0.  2  sec/cm  n  and  the  ordinate  for  Zlog 
signal  is  one  decade  of  irradiance  per  cm.  The  linear  baseline  is  one 

cm  up  from  the  bottom  of  each  frame,  and  the  log  baseline  is  indefinite 
ana  not  pertinent. 


15.  Same  as  Figure  14,  for  10.  6  jim. 
Figure  8. 


The  turbulence  levels  are  those  of 


16.  Cumulative  probability  distributions  for  log  irradiance. 
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Curve 

Wavelength 

C  2 

n 

A 

4880# 

4.  5  x  10 

B 

4880^ 

3.4  x  10“12 

C 

10.  6  [lm 

4.2  x  10'15 

D 

1 0.  6  am 

3.2  x  10'12 

For  curve  B,  the  abscissa  should  be  multiplied  bv  3  3  The  h  i, 
curve  represents  a  Rayleigh  amplitude  distributton. 


17.  Illustration  of  the  application  of  reciprocity. 


(a)  Target  illumination  system 

(b)  Reciprocal  (conceptual)  heterodyne  receiver  or 
imaging  system 

19c  Mean  target  irradiance,  normalized  by  that  for  a  verv  br 

as  a  function  of  aperture  (after  Ref.  25).  7  **  aperture» 

20‘  var„l0  vCs.diaap8errre°f  ^radiance 

various  regions  of  (b/rj  are  indiCIII.  WMCh  pred°m“ate  f- 

vs.  aperture,  for  beam  winder  awl.^  ****  "°rmalized  irradia"ce  varianc 

22‘  flc^tefnlmiueri  ^  •*  °f  * 

distance  z1  fmm  n  ,  ..  L’  The  scatterer  is  located  at  a 

resect  te  theTrig^att;  "T  ^  a*  ^  Sm*n  “*ta  ^ 

”*  collimated* transmitter  **  a  ^r-fieid. 

transmitter.  throughTsmaU  Ingle  “  dU*“”  *“  «“ 

24.  Illustrating  a  ray  which  is  scattered  out  of  the  focal  snnf  in  r- 

which  can  interfere  with  unperturbed  rays  in  the  far  field  of  R.8^6  ' 

25*  unnStrt^g/Cattering  int°  R  from  an  ^"Perturbed  ray  out  of  R  The 
unperturbed  ray  arises  due  to  aberrations  or  defocus ine  The  fh  !  , 

ray  ts  at  an  angle  v  with  respect  to  tee  optical^Tn'd  te^scattetr 
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shown  is  at  midpath  with  a  scattering  angle  (3. 

2h.  Plane  wave  propagation  through  a  slab  of  uniform  turbulence,  representing 
an  idealization  of  turbulence  intermittency. 

27.  Illustration  of  random  process  made  up  of  the  product  of  fast  and  slow 
subprocesses. 

28.  Power  spectra  of  the  processes  of  Figure  27. 

29.  Square  of  the  random  process  (x)  of  Figure  27. 

30.  Power  spectrum  of  the  process  of  Figure  29. 
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(Aperture  Diffraction) 


FIGURE  18 


FIGURE  20  log  b 


Region  of  Scintillation 


FIGURE  24 


FIGURE  25 


